Spin relaxation in quantum Hall systems 
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We study the spin relaxation in an interacting two-dimensional electron gas in a strong magnetic 
field for the case that the electron density is close to filling just one Landau sub-level of one spin 
projection, i.e., for filling factor f ~ 1. Assuming the relaxation to be caused by scattering with 
phonons, we derive the kinetic equations for the electron's spin-density which replace the Bloch 
equations in our case. These equations are non-linear and their solution depends crucially on the 
filling factor and on the temperature T of the phonon bath. In the limit of T = and = 1, the 
solution relaxes asymptotically with a power law inversely proportional to time, instead of following 
the conventional exponential behavior. 
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Experiments in two-dimensional (2-d) electronic struc- 
tures in a strong magnetic field reveal a rich variety of 
features. Recently, many experiments have focussed on 
properties related to the electrons' spins Optical 
measurements of the magnetization and transport mea- 
surements of the activation energy are interpreted with 
the concept of Skyrmion-quasiparticles |^Jl^ . These can 
be described as extended spin-textures containing a num- 
ber of spins which are flipped with respect to the prefer- 
ential direction set by the magnetic field. Their number 
is determined by the competition between Zeeman en- 
ergy and Coulomb interaction of the electrons. While the 
equilibrium of systems containing Skyrmions has been 
studied in various ways - both experimentally and theo- 
retically, the Skyrmion's dynamical properties are much 
less understood. There are indirect measurements of elec- 
tron spin resonance which use the magnetoresisitivity to 
determine the spin-splitting of the Landau levels 11-1^] 
and which also yield a line width |p^ . 
The first and fundamental question about the dynamics 
regards the formation of the Skyrmion. A Skyrmion anti- 
Skyrmion pair has only half the energy of a particle-hole 
excitation, a spin exciton [9|jl0|]; nevertheless, it needs 
for its formation a mechanism which flips the electron's 
spins, since a Skyrmion may contain a large number of 
overturned spins. The mechanism must be provided from 
outside the bare electron system of Coulomb interaction 
and Zeeman energy, because the spin parallel to the mag- 
netic field is a constant of motion. In this work, we wish 
to analyze the spin relaxation, originating from such a 
mechanism, specifically for a system of 2-d interacting 
electrons in a strong magnetic field. As it will be shown 
below, the result is intriguing and non-trivial for this 
system, already in the case of a uniform spin-density, 
and quite different from the well-known Bloch equations 
0. Thus, we will demonstrate our new approach study- 
ing the uniform case; this has to be seen as the first step 
towards a solution of the more general problem of the 
formation of the Skyrmion. 

We consider in the following interacting electrons in two 
dimensions moving in a strong magnetic field under the 



following conditions: the electron density should be close 
to that of one filled spin-split Landau sub-level, i.e. the 
filling factor v ^ 1, and the Zeeman energy should be 
much less than the Coulomb energy. Both conditions are 
relevant for the Skyrmion case. The latter is realized 
in the experiment by a variation of the ij-factor with 
pressure In three-dimensional (3-d) semiconductors, 
various mechanisms of spin relaxation have been studied 
experimentally (e.g., in optical measurements) and theo- 
retically. Among them are scattering by holes due to the 
exchange energy, and scattering by phonons and impu- 
rities which becomes possible because of the spin-orbit 
interaction; for a comparison of these different mecha- 
nisms, see Ref. | p5| . In a system in strong magnetic field, 
the Zeeman energy has to be taken up by the scatterer 
in an elementary spin flip. Therefore, we study here a 
general model of inelastic scattering by phonons. In de- 
riving the kinetic equation which governs the relaxation 
of the spin-density due to the scattering by phonons, 
the following three assumptions are made: (i) the 2-d 
electrons are described by a non-equilibrium distribution 
which yields uniform expectation values for density and 
spin-density, (ii) the phonons are in thermal equilibrium 
unperturbed by the electrons, (iii) the electron-phonon 
coupling is so weak that it can be treated in perturba- 
tion theory. The appropriate technique for performing 
perturbation theory in non-equilibrium systems was in- 
troduced by Keldysh |16 1^. It is shown in Ref. |l^] 
how to derive the kinetic equation in the quasi-classical 
regime. The condition for the applicability of the quasi- 
classical approximation in the present context of the spin 
relaxation is given below. 

How can phonons flip the spin ? Consider the standard 
expression of the spin-orbit energy. 
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where ^cr and p are spin and momentum of the elec- 
tron. The electric field E has its origin in a piezoelectric 
distortion of the lattice and it is linear in phonon opera- 
tors. Thus, the electrons are scattered and their spin is 



1 



flipped by the phonons of the 3-d lattice which contains 
the 2-d electron system at z = 0. The orbital part of the 
electron's single particle states is confined to the lowest 
Landau level. We denote the projected spin-density of 
the electrons at (2-d) wave vector q with S(q) Then, 
we get for the electron-phonon part of the Hamiltonian 

6q J creates a phonon with (3-d) wave vector Q ~ (q, Q') 
and polarization s. The coupling constant As(Q) can 
be derived from (|l|). It containcs the common electronic 
wave function in the z-direction. With a relation between 
electric field and phonons as in the model of Refs. p9| , p0| , 
As(Q) is quadratic in momentum. The considerations 
below are valid for arbitrary As(Q). 
Having specified the model, we will now study first the 
kinetic equation for a strongly simplified case in which 
only two non-interacting electronic states are kept. In 
a second step, we will then come back to the full model 
of interacting electrons in two Landau sub-levels coupled 
to a phonon bath. 

Two non-interacting states — The Hamiltonian is given 

by 

Htoy= - A(c|c| -cjcJ+^Wq ffo]j6q+ij 

+ Y.\{clc^b^ + c\c^bl) , (3) 
q 

where c| (6]j) are electron (phonon) creation operators. 
Here, scattering from the t state to the J, state occurs un- 
der absorption of a phonon with an arbitrary wavevector 
q, and the same is true for the reverse process. The 
kinetic equation for the time-dependent average occu- 
pation numbers n|^|(t) can be derived with the aid of 
Fermi's golden rule (here, time arguments are omitted): 

dtn^ = -dtn^ - - [(1 - n^hiil + N)-n^il- n^)N] . 

T 

(4) 

Eq. (^) is valid in the quasi-stationary limit A >> 1/r, 
and the relaxation time r is given by 

^ = J2xl2n6{2A-oo^). (5) 
q 

iV = [e^^'^ ~ 1] ^^ is the occupation number of the only 
phonon {ujq — 2A) which can be effective under energy 
conservation. The density v = ni{t) + ni{t) is fixed; thus 
it is clear that Eq. (|^) has two stationary solutions for 
n|. Only one of these is physical, the other is unphysi- 
cal (n| > 1). The kinetic equation (0) is best discussed 



in terms of the depolarization, the deviation from the 
physical solution, d{t) = n|(cx)) — rit(t) : 

r dt S{t) = -S{t) [e + S{t)] . (6) 

The crucial parameter e is the difference between the two 
stationary solutions, e = [4iV(l + N) + {ly ~ 1)^] It 
depends strongly on the temperature and the total den- 
sity. For zero temperature 1//3 = 0, e = 1|. For = 1, 
on the other hand, e ^ 2e~^^ at low temperatures. Now 
it is obvious that, and under which conditions, lineariza- 
tion in the kinetic equation (^) can fail: Eq. (^ shows 
two regimes depending on the relative size of e and 5{t). 
If 6{t) is smaller than e, then one can linearize Eq. (^ in 
5{t), and 5{t) decays exponentially. If 5{t) is larger than 
e, then Eq. (^ becomes quadratic in 5{t), and 5{t) decays 
with \/t until it becomes so small that it crosses over into 
the first regime. This behavior is pictured in Fig. 1 for 
5{t = 0) = 0.25. The solid line shows the depolariza- 
tion for e = 0.05. Initially, the curve follows that of the 
quadratic equation (e = 0, dashed line), until it crosses 
over to exponential behavior. This limiting exponential 
behavior is indicated as a dotted line. We conclude that 
even the solution of the oversimplified model Htoy shows 
a non-trivial behavior, far from being exponential per 
se. It displays a strong dependence on temperature and 
density. The asymptotic exponential decay of 5{t) can 
become arbitrarily slow (e 0) as the temperature ap- 
proaches and — > 1; and in the limit, 5{t) oc \/t. 
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Fig. 1: Depolarization &{t) as a function oi t/r for an 
initial value of (5(0) — 0.25; solid line: e = 0.05, dashed 
line: e = 0. 

Electron-electron interaction — We now turn back to 
the case of interacting electrons in the two lowest Lan- 
dau sub-levels split by the Zeeman energy. The electrons 
are scattered by phonons, see He-ph, Eq. (|^). The elec- 
tronic part He of the total Hamiltonian He+Hph+He-ph 
is then 

He ^ - A ^(cjpC^p - cjpC^p) + Hcoul ■ (7) 

p 

Here, cj^ (cjp) creates a Landau state with linear mo- 
mentum p and spin ] {[). The most important ingredi- 
ent, which was missing in Htoy, is the Coulomb interac- 
tion Hcoui between the electrons, cf. pll. As stated in 
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assumption (i) above, the expectation values of density 
and spin-density of the electrons are presumed to be uni- 
form. The density is conserved and its value is given by 
the filling factor ~ 1. We wish to study the uniform 
spin-density, S(t) =< S(q = 0) >, 



-/V$ denotes the number of states in a Landau sub-level; 
if all I states are populated and all J, states are empty, 



The derivation of the kinetic equation 



for S(t) can be performed with the method of Keldysh 
| p6[ . The kinetic equation contains a precession term, 
which is already present without any phonon scattering, 
and the main term resulting from the collisions with the 
phonons. In leading (second) order perturbation theory 
in the clectron-phonon coupling, the collision integral be- 
comes 



Q j.kd 



[C'=' (q; i, t')V=' {Cl-t-t')-C'^ (q; t' , t)V'^ (Q; t' - t)] . 

(9) 

Only times t' < t contribute due to causality. The func- 
tion C is the dynamical spin-spin correlation function, 

C^-'(q;i,t') = ^ < S\q,t) S\-q,t') > . (10) 

The spin-density couples in H^-ph to the phonons via 



*(Q) 



(11) 



thus, the phonon expectation value I? in (|9[) is given by 

I?J'(Q;<-i') =< ^^(Q,t) ^'(-Q,t') > ■ (12) 

D is easily calculated with the equilibrium Hamiltonian 
Hph of the phonons. e*-''^ is the antisymmetric tensor. 
Its origin is the time derivative in (^) which leads to a 
commutator of S with the spin-density in TJe-p/i- 
It is apparent now that collective modes are responsi- 
ble for the relaxation process. The question is whether 
one can express the two-particle (four-fermion) correla- 
tion function of the interacting system, C, again by S(t) 
and thus derive a closed equation for S(i). Fortunately, 
that turns out to be possible in our case: The time de- 
pendence of C is determined by particle-hole excitations 
(e.g., S+ = + iSy - c\ci). For H^, the particle- 
hole excitations above the ground state are spin-excitons 
and their dispersion E^x (q) is rigorously known |^,^ at 
^ = 1, 
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(13) 



K is the dielectric constant, Ib is the magnetic length, and 
Io{x) denotes the Bessel function. We now approximate, 
in the present case of a non-equilibrium state and z/ ~ 1, 
the time dependence of C with the above dispersion of 
the spin-exciton, Eex{q)- Thus we neglect (i) the differ- 
ence between the real energy of a single spin-exciton in 
a non-equilibrium state and E^x (q) , and (ii) interactions 
between the spin-excitations at v ^ 1. This yields, e.g.. 



C'^'+(q;t,t') ~ e*^"(«)(*-*') C''+{ci;t,t) 



(14) 



The remaining equal-time correlations are then calcu- 
lated in the Hartree-Fock approximation. This is the 
weak-scattering approximation in which the effect of the 
phonons is included in the lowest order in the state, but 
neglected in the time dependence 0]. Since the char- 
acteristic energy for the time dependence is the Zeeman 
energy A, that corresponds to the condition A >> l/r, 
where the scattering time r is defined below in Eq. (|l^) . 
The two approximations above regarding the exciton en- 
ergy demand that the electronic temperature must be 
small compared to the Zeeman energy A and also, the 
filling factor must be such that — 1| << 1. Now, col- 
lecting all terms, we get the resulting kinetic equations 
for the components of S{t) {S~^ = -\-iSy\ details of the 
calculation arc deferred to a forthcoming publication): 



s'-{t) 



'-NS'{t) 



(15) 



{dt+i2l\) S+it) 



S+it) 



The relaxation time r and the average phonon number 
N are defined by the two components of the following 



equation (here, N{u}) 



l]-i) 



l/r 
N/t 



j = x,y 



1 

(16) 



It is quite instructive to formulate again the kinetic equa- 
tion in terms of the depolarization vector S{t). We use 
(5^(t) = S^{oo) — S^{t) and we split the precession term 
off S+{t) by redefining 5+{t) = cxp(i2At)S'+(t). Then, 
the result is 



dtS+it) 



(17) 



Here, the parameter e corresponding to the one used 
above in the case of the two-state model is defined as 
e = [4iV(l + N) + {iy~ 1)^] Obviously, Eqs. (|l^ are 
quite different from the standard Bloch equations p4|. 
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They are non-linear and it is seen that the relaxation 
of the transverse component 6^ (t) depends on the other 
component. Since the first of the two kinetic equations 
( p7| ) is identical to (||), the same discussion applies. Ob- 
serving TdtHS''it)/6+{t)] = -e/2, Eqs. @ can be ex- 
plicitely solved with the result 



(5^(0) 



- p-Et/(2r) 



-ttjT 



(5+(0) 1 + (5^(0)(1 -e-<=*/^)/e 



(18) 



For any finite value of the parameter e, the leading 
asymptotic behavior yields a ratio T2 = 2Ti of the re- 
laxation times. On the other hand, for zero temperature 
and = 1, all normalized components of the depolar- 
ization follow the same function which is a power law 
oc t-^. 

The relaxation time r is determined by those phonons, 
whose energy and in-plane momentum match the energy 
and the momentum of the spin-excitons, see Eq. (16). 
Using the model of D'yakonov and Perel' (c.f. Ref. 
1^,^), we calculate the electron phonon coupling pa- 
rameter A^(Q) and perform the summations and inte- 
grations in (16) with the following result in the limit of 
A = 



1 _l [2 /i;n2 h 
r ^ 8 V ^ Vs/ p/| 



(19) 



Here, v parametrizes the electron phonon coupling, s 
is the phonon velocity, p the 3-d density, and Wi.i{x) 
the Whittaker function. The parameter xo is given by 
xq = {hs/{eclB))'^ , where ec is the Coulomb energy. An 
estimate with characteristic values for electrons in GaAs 
at a magnetic field of 10 T yields r ~ 10~^*'sec. This 
order of magnitude is much smaller than the estimate in 
Ref. po| and in agreement with experimental observation 

We have studied the spin relaxation in a quantum Hall 
system for filling factors ~ 1. The resulting kinetic 
equation ( [l5| ) is non-linear and quite different from a 
conventional Bloch equation. This difference is most pro- 
nounced for v —^ 1 and T — > 0; there, a linearization of 
the kinetic equation fails, since the linear term vanishes 
at = 1, T = 0. The relevant time r (|l|,|l^ in the ki- 
netic equation is determined dominantly by the Coulomb 
interaction between the electrons. — 
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